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Abstract The conditions for the existence and stability of cosmological power-law scaling solutions are es- 
tablished when the Einstein-Hilbert action is modified by the inclusion of a function of the Gauss-Bonnet 
curvature invariant. The general form of the action that leads to such solutions is determined for the case 
where the universe is sourced by a barotropic perfect fluid. It is shown by employing an equivalence between 
the Gauss-Bonnet action and a scalar-tensor theory of gravity that the cosmological field equations can be 
written as a plane autonomous system. It is found that stable scaling solutions exist when the parameters of 
the model take appropriate values. 
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1 Introduction 

In recent years, there has been considerable interest in the possibility that Einstein's theory of general relativ- 
ity may become modified in high-curvature regimes and over large distance scales. This possibility has been 
motivated by a wealth of high redshift observations, which indicate that the universe is presently undergoing 
a phase of accelerated expansion lTll[^f^mi5ll^mf8ll^rT0l[TTirrarT3irrarT^ . (For a review, see, e.g. ED). In 
many such studies, the Einstein-Hilbert action is modified by the introduction of terms involving higher- 
order curvature invariants. An important quadratic combination of such invariants, which is motivated by 
string theory, is given by the Gauss-Bonnet (GB) invariant 

^= R 2_ 4R ^v R ^ +R nv P r R ^ (1) 

In four dimensions, the GB term is a topological invariant and introducing a term proportional to into the 
Einstein-Hilbert action does not modify the dynamics. Recently, however, the cosmology of models based 
on a class of generalised theories with an action of the form 

(2) 
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has been considered, where ) is a differentiable function of 5f and 5 m represents the matter action |[T7l 
[Ill[l9l|20l|2ll[22ll23l[2ll[25l|2^ . 

The purpose of the present paper is to investigate the existence and stability of cosmological power-law 
scaling solutions derived from theories of the type © in the presence of a perfect fluid matter source. Scaling 
(attractor) solutions play an important role in cosmology, since they enable the asymptotic behaviour and 
stability of a particular cosmological background to be determined. Moreover, they provide a framework for 
establishing the behaviour of more general cosmological solutions H321l331l34ll351l361l37l . 

The structure of the paper is as follows. We begin in Section 2 by summarizing the derivation of the 
cosmological field equations by employing an equivalence between the action © and a corresponding action 
involving a self-interacting scalar field that is non-minimally coupled to gravity. We focus on the spatially flat 
and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) universe and proceed in Section 3 to identify 
the most general form for the function /(5f ) that results in power-law (scaling) solutions when the matter 
source is a barotropic fluid with a constant equation of state parameter. Specifically, we find that scaling 
solutions may arise when / = ±2V CttS, where a is an arbitrary constant. We then show that for this form 
of the action, the field equations can be expressed as a plane autonomous system. This allows us to employ 
dynamical systems theory to investigate the stability of the vacuum and non-vacuum solutions and this is 
done in Sections 4 and 5, respectively. We find that scaling solutions, corresponding either to a stable node or 
a stable spiral node, can arise when the equation of state of the fluid and the parameter, a, satisfy appropriate 
conditions. We conclude with a discussion in Section 6. Units are chosen such that SnG = c = 1. 



2 Cosmological Field Equations 

Action © may be expressed in an alternative form by introducing two auxiliary scalar fields % and £ such 
that ifTTir 



s = | (~ + C(^ -*)+/(*)) +s» 



(3) 



Varying Eq. © with respect to £ yields the constraint % = thereby reproducing action ©. On the other 
hand, varying action © with respect to % implies that £ =F(x)> whereF(^) = d f{%) / d%, and substituting 
this condition back into Eq. © leads to 

5 = Jd 4 x^(^+F( % )(y-x)+f(x)^+S m . (4) 

It follows, therefore, that the action © is equivalent to the action 11171,138] 

S = J rfV=g(~-V (0)- h{4>)<^+S m (5) 

where the scalar field, 0, is defined implicitly by 

h ft) = _F(^) (6) 
for some function h(<j>) and has an effective self-interaction potential 

Vf» = SfF(Sf)-/(Sf), (7) 

where F = df /dW. Eq. © may be interpreted as an effective 'scalar-tensor' theory, where the scalar field 
has a vanishing kinetic term. 

To study cosmological models based on action ©, one may proceed directly by varying the action to 
derive the field equations or, indirectly, by varying the equivalent action ©. We employ the latter approach 
in the present work in view of its potential simplicity. The field equations in this case take the form 

R HV_}_ Rg HV = T £v + T gv (8) 
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where T% v is the energy-momentum tensor of the matter fields and T^ v denotes the effective energy- 
momentum tensor resulting from the scalar field, 0, and the GB term. Since the GB term is a topological 
invariant in four dimensions, the standard field equations of GR are recovered when h((j>) = constant. Con- 
sequently, only terms involving derivatives of h((j>) arise in the energy-momentum tensor, which is given by 
ED 

= - ^ v y(0)-2[V^V v / J (0)]/? + 2^ v [V 2 / J (0)]/? + 4[V p V^/ J (0)]7? v " 
+ 4[V p V v fc(0)] J RW - 4[V 2 /*(0)]/^ v - 4^ v [V p V T /<0)]/?" T 

+ 4[V p V T /j(0pw VT . (9) 

Finally, the equation of motion for the scalar field takes the form 

V (0) +h t t(<t)9 = 0, (10) 

where a comma denotes differentiation with respect to 0. 

Our aim is to study the dynamics of the isotropic and spatially flat FLRW universe sourced by a perfect 
barotropic fluid with an equation of state parameter, w m = p m / p m , where p m and p m denote the pressure and 
energy density of the fluid, respectively. For this spacetime, the GB invariant is given by <3 = 24H 2 (H +H 2 ), 
where H = d/a defines the Hubble parameter, a represents the scale factor of the universe and a dot denotes 
differentiation with respect to cosmic time. The Friedmann and Raychaudhuri equations derived from Eqs. 
©-© for this background are then given by 11911421 



3H 2 = 7(0) + 24// 3 A + p m , (11) 

^ + 3^ H 2 = 7(0) + m 2 h + 16// 3 A (l + - p m , (12) 

respectively, and the scalar field equation ( fTOl reduces to 

V^+24h A H 2 (H + H 2 )=0. (13) 

It proves convenient to interpret the GB gravitational terms on the right-hand side of the Friedmann 

equation (fTTT ) as an effective energy density, such that p$ = Ty + 7(0), where % = 24hH 3 plays the role of 
a kinetic energy. It is then natural to introduce the dimensionless variables 

7(0) 

= ~o77T> yi = SHh, (14) 



3H 2 



and the fractional energy densities 



an = ^ = i-yi-yi, (15) 

G&^yi+yi- (16) 

The background field equations (ITTb-(fT3T> can then be expressed in terms of these variables such that 

d ^ = 2ey l -(\-e)y 2 , (17) 

^ = -2s + 3(\ - yi ) - (2- s)y 2 + 3w m n m , (18) 



where £ = 



-H/H 2 and N = In a. 
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3 Cosmological Scaling Solutions 

We wish to identify the class of GB theories that admit scaling solutions such that each of the terms in the 
Friedmann equation (TTTb scales at the same rate, H 2 

x Pm x ^(0) °^ 1143 H « These conditions result in a 
power-law solution to Eqs. (fTTI) -(fT3l) of the form a°^t l l e , where e = constant. For such a scaling solution, 
it follows from Eq. (fT3l) that 

V 4 = -^V 2 h^ (19) 

when £ / 1 , where a is a finite constant. Integrating Eq. (fT9l) then implies that 

oc 

h=y+P, (20) 

where /3 is an arbitrary integration constant. 

Relating the functions V(<p) and h(<p) in this way is equivalent to specifying the form of the GB function, 
via the definition given in Eq. |7]). Indeed, substituting Eq. (1201) into Eq. (0) results in the first-order, 
non-linear differential equation 

,<f-A(<f + S\ = -*. p „ 



d<S J )\d ( S 

Eq. (1211) is an example of Clairaut's equation l44l and may be solved in full generality by differentiating 
with respect to ( S: 

d 2 f \ ( df n \ 2 a 

+ P ) ~& =0 ' (22) 

Eq. (122)) is trivially solved by ) = CCo + a0, where a,- are constants. However, this simply corresponds 
to the introduction of a cosmological constant in the action (0 and is not physically interesting to the 
present discussion. (Recall that a contribution of the form / <=c <~g is also uninteresting since the GB term is a 
topological invariant). On the other hand, a singular solution to Eq. (1211 ) with no arbitrary constants can be 
found by setting the square bracketed term in Eq. (1221 to zero and substituting the result into Eq. (l2Tb . We 

find that 

/(Sf ) = ±2Va&, (23) 

where we have specified j3 = without loss of generality. Moreover, requiring the action (fSJ) to be real 
implies that a5f > 0. 

Eqs. (20) and (23) represent the necessary and sufficient conditions for the existence of power-law scaling 
solutions, where e = constant. More general solutions to the field equations, where e is time-dependent, exist 
for this model. If the cosmological behaviour of the model (1231 is to be determined, the coupled differential 
equations (fT71)-(fT8T) must close. This implies that the parameter e must be expressible as a function of yi and 
y>2 only. When Eq. (1251 is satisfied, we find that 

— l-E* P4) 
Hence, substituting Eq. (1241 ) into Eqs. (fT7T)-(fT8l) yields the plane autonomous system: 

j| = 2(j 2 - 1 ) - ^Ljfo + ^ j? + 3 ( 1 + w m ) ( 1 - yi - w ) . (26) 

Before concluding this section, it should be remarked that the equivalence between actions © and © 
does not apply for the special case e = 1 (yi =0), corresponding to the coasting solution, a <=c t. In this case, 
integration of Eq. (fT3l would yield V((f)) = Vo = constant and the solution to Eq. |7]) would then be given 
by /(Sf ) = —Vo + for some constant y. This disparity can be traced to the singular nature of the coasting 
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solution for the model (|23l . Specifically, the Friedmann equation derived directly from action Q for this 
model is given by 

, , — H 2 (2H 3 -H) 
3// 2 = T Vto ^ +H2? ,, +P m (27) 

and the term originating from the GB contribution is ill-defined when e = 1 (y\ = 0). Consequently, we do 
not consider this solution in the phase plane analyses of the following sections. 



4 Vacuum solutions 

In this Section, we consider vacuum solutions where Q m = and y\ =l—y 2 . The pair of equations d25T)-(l26l> 
then reduces to the one-dimensional system 



There exist two power-law solutions when y\ ^ 0: 

yi = _I±IV9 + I92a, (29) 
2 o 

which we denote as "V , respectively. The reality of the fixed points requires that a > —9/192. The power 
of the expansion can be expressed in terms of the effective equation of state parameter 

w e// =-l + |e (30) 

such that a(t) oc f 2 /[ 3 ( 1 + H V/)l. it is determined by the value of the GB coupling parameter, a, and substituting 
Eqs. (|24]> and ([29]> into Eq. ([30]) implies that 



W eff 



24a 



-40a - 3 ± V9 + 192a 



(31) 



where the +/— corresponds to the points Y ± , respectively. This dependency of the effective equation of 
state on the GB parameter is illustrated in Fig. [TJ The solution Y + corresponds to an inflationary cosmology 
when a > and the exponential, de Sitter solution arises when a = 3/8. The solution "V~ is in a super- 
inflationary regime (w e // < —1) for a > 0. When a < 0, the effective equation of state corresponds to that 
of an ultra-stiff fluid (w e ff > 1). Our results are in line with the recent conclusions of Ref. [45], where a 
study of the late-time cosmology based on the model f(&) ^ — @" was made with the field equations derived 
directly from action (f2]). 

The eigenvalues associated with the equilibrium points "V are given by 

M ± = -4 - tJ- ± -j- V9 + 192a. (32) 
16a 16a 

The solution 7 /+ is stable for a > —9/192. The solution f~ is a stable point when a > and unstable for 
-9/192 < a <0. 
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Fig. 1 Illustrating the effective equation of state w e ff for the vacuum solutions */ /+ and Y~ . Requiring that the fixed points are real yields the 
condition a > —9/192. The left-hand panel corresponds to */ /+ , which shows that Y + is an accelerating solution for a > and corresponds to the 
de Sitter solution if a = 3/8. The middle panel corresponds to Y~ when a < and in this regime w e ff > 1. The right-hand panel corresponds to 
"V~ when a > and in this regime w e ff < — 1. 



5 Non-vacuum solutions 



In this Section, we study the background dynamics of models based on GB theories of the type (1231 ) in 
the presence of a perfect fluid. The vacuum solutions remain as equilibrium points of the autonomous 
system (|25T>-(|261): 

{yuyi)= f-i±iv9 + 192a, \t^V9 + 192a). (33) 



In addition, there exist two scaling solutions, where Q m and £hg are constants: 

12a(l+w 



&2n 



2 v / -3a(l + 3w„ 



IT 



2 v / -3a(l + 3w M 



2 v /-3a(l + 3w„ 



^-3a(l + 3w m 
12a(l +w m ) 
v / -3a(l + 3w m )' 
12a(l +w m ) 




^-3a(l + 3 



and w e ff = w m . We denote these solutions by 5? . 

The eigenvalues associated with the equilibrium points are given by 



l 



32a 
1 

32a 



48a(3 + w m ) + 9 =f 3V9 + 192a +X{ 



48a(3 + w m ) + 9 =f 3 V9 + 192a -X[ 



(34) 
(35) 
(36) 



(37) 



(38) 



32a 



256a 2 (l + 3w m ) 2 + 288a(l + w m ) + 18 =f 32a(l + 3w m )V9 + 192a t 6V9+192a 



1/2 



(39) 



The stability of these vacuum solutions is altered when a matter source is introduced into the system and 
depends on both the GB parameter, a, and the perfect fluid equation of state, w m . This dependency is 
illustrated in Fig. [2 The solid lines represent the regions where the nature of the equilibrium points changes 
as the parameter values are altered. The stability of y~ is determined by the sign of the GB parameter, a. 
On the boundary distinguishing the nature of the fixed point ~f + , one of the eigenvalues ju^ 2 vanishes. To 
analyse the stability of the equilibrium point for these particular choices of parameter values would require 
a second-order analysis, which is beyond the scope of the present work. 
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Fig. 2 Illustrating the nature of the equilibrium points Y + (left-hand panel) and (right-hand panel) in the parameter space spanned by (w,„, a). 
Both fixed points are real if a > —9/192. On the boundary (denoted by the solid line) that distinguishes the stability of the fixed point 'Y^ , one of 
the eigenvalues jU[t, vanishes. This is indicated in the figure by a change in colour. The dotted line in the left-hand panel represents the invariant 
sub-manifold yi = 0. In the case of the point S~ , neither of the eigenvalues vanishes in any region of the (w m ,a) plane. 



The eigenvalues associated with the scaling equilibrium points y ± are given by 
3 

4« 

3 



-(w„ 



1) + ^ 



A*2 



1 

4a 



1 



±8a(l + 3w m ) v / -3a(l + 3 



w„ 



a 2 (135w,^ + 306w m 



+ 71) 



1/2 



(40) 
(41) 
(42) 



The stability of these fixed points is illustrated in Fig. [3] The points are real in the region of parameter space, 
ce(l + 3w m ) < 0. Furthermore, they are only physically meaningful if Q m = 1 — y\ — y2 > 0. This results in 
a further restriction in the (w m , a) plane after substitution of Eq. (l35l) . 

The top two panels of Fig. [3] correspond to the scaling solution S" + where y\ > and the bottom two 
panels correspond to y~ where y\ < 0. The point S" + is either a stable node or a stable spiral. The point 
is always a saddle. On the curve one of the eigenvalues of vanishes. 

To illustrate the scaling dynamics, let us consider the specific case where (a,w m ) = (0.05,-0.6). At 
this location in parameter space, there exist two equilibrium pointg[J the saddle point Y + and the stable 
node y + . The basin of attraction for S p+ is shown in Fig. [4] As a second example, we consider the case 
(cc,w m ) = (—0.005,-0.05), where there exist four equilibrium points: an unstable vacuum solution a 
saddle point a stable "V + and a stable spiral The spiral nature of the point ^ + is illustrated in the 
phase portrait of Fig. O where the initial conditions were specified to be Q m = = 0.5. 



6 Conclusion 

In this paper we have investigated the existence and stability of cosmological power-law scaling solutions 
sourced by a barotropic fluid when an appropriate function of the Gauss-Bonnett topological invariant is 
introduced into the Einstein-Hilbert action. It was found that the general class of such theories that ad- 
mit power-law solutions is given by Eq. d23l) , i.e., f(&) = ±2V a'S for some constant coefficient, a. By 
exploiting an equivalence between generalized Gauss-Bonnet gravitational theories and a corresponding 
higher-order, scalar-tensor theory, it was further shown that the Friedmann equations for this class of model 
can be written in the form of a two-dimensional dynamical system. The stability of the equilibrium points for 

1 Note that the point also exists but this occurs in the region y\ < 0. Stable scaling solutions arise only for y\ > and, since y\ = is a 
separatrix, a trajectory beginning in the region yi < will not be able to reach . We therefore choose the initial conditions in Fig.f4]such that 
y\ > 0. This is equivalent to choosing the negative sign in Eq. (23). 
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Fig. 3 Illustrating the stability of the scaling equilibrium points ,9^ in the parameter space spanned by (w m , a). The region of parameter space 
is restricted by the requirement that the equilibrium points are real, a(l +3w m ) < 0, and also correspond to physically realistic solutions where 
&m > 0- The shaded areas depict the regions of parameter space where the solutions are unphysical. These restrictions imply that the analysis can 
be separated into regions where a > (left-hand panels) and a < (right-hand panels). The regions of parameter space where the fixed points 
correspond to either a saddle point or a stable/spiral node are identified. On the line Q m = 0, the eigenvalue flf = (for the scaling point 
when a > 0. Conversely, ;U^~ = (for the scaling point J? - ) when £2„, = and a < 0. 




both vacuum and non-vacuum models was established. In the former case, the GB parameter, a, determines 
the effective equation of state parameter. For non- vacuum solutions, the nature of the critical points depends 
on both a and the fluid equation of state parameter, w m . The regions of parameter space (a,w m ) that admit 
stable non-vacuum scaling solutions were identified. 

The models we have investigated do not admit a transition from a decelerating to an accelerating phase 
of cosmic expansion. However, our aim in this paper has been to focus on power-law solutions rather than 
develop a phenomenological model of generalized Gauss-Bonnet gravity as a candidate for dark energy. 
Power-law solutions are of interest since they can be regarded as approximations to more realistic models. In 
particular, phenomenological models could be constructed where the parameter a is given by some function 
of (or equivalently the scalar field 0), such that a is slowly varying for much of the history of the universe, 
but at some epoch undergoes a change in sign. In principle, this could cause the universe to enter a phase 
of accelerated expansion. It would be interesting to develop specific models of this type, along the lines 
outlined in Ref. W\ . 
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Fig. 4 Illustrating the dynamics of the model j23t for the particular case where (a,w,„) = (0.05,-0.6). The left-hand panel depicts the phase 
space, where the straight line y\ = 1 — yi corresponds to the vacuum solution £2 m = 0. The red dot represents the scaling fixed point ,y + . For the 
range of initial conditions chosen, all non- vacuum, physically acceptable solutions are attracted to .9' + . The right-hand panel depicts the evolution 
of the fractional energy densities of the perfect fluid, Q. m , and the GB contribution, Cluj, for the initial conditions i2„, = £2cf = 0.5. It is seen that the 
fractional densities asymptote to constant values at late times, thus indicating that the solution is scaling. 




Fig. 5 Illustrating the dynamics of the model {23) for the particular case where (a,w m ) = (—0.005, —0.05). The left-hand panel depicts the phase 
space for this scenario, whereas the right-hand panel depicts the evolution of the fractional energy densities Q„, and 12<y. The initial conditions were 
chosen such that Cl m = Qcg = 0.5. At late times, the fractional energy densities of the fluid and GB contribution tend to constant values. 
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